Natural variation in neuraminidase activity influences the evolutionary potential of the seasonal H1N1 lineage hemagglutinin.

The antigenic evolution of the influenza A virus hemagglutinin (HA) gene poses a major challenge for the development of vaccines capable of eliciting long-term protection. Prior efforts to understand the mechanisms that govern viral antigenic evolution mainly focus on HA in isolation, ignoring the fact that HA must act in concert with the viral neuraminidase (NA) during replication and spread. Numerous studies have demonstrated that the degree to which the receptor binding avidity of HA and receptor cleaving activity of NA are balanced with each other influences overall viral fitness. We recently showed that changes in NA activity can significantly alter the mutational fitness landscape of HA in the context of a lab-adapted virus strain. Here, we test whether natural variation in relative NA activity can influence the evolutionary potential of HA in the context of the seasonal H1N1 lineage (pdmH1N1) that has circulated in humans since the 2009 pandemic. We observed substantial variation in the relative activities of NA proteins encoded by a panel of H1N1 vaccine strains isolated between 2009 and 2019. We comprehensively assessed the effect of NA background on the HA mutational fitness landscape in the circulating pdmH1N1 lineage using deep mutational scanning and observed pronounced epistasis between NA and residues in or near the receptor binding site of HA. To determine whether NA variation could influence the antigenic evolution of HA, we performed neutralizing antibody selection experiments using a panel of monoclonal antibodies targeting different HA epitopes. We found that the specific antibody escape profiles of HA were highly contingent upon NA background. Overall, our results indicate that natural variation in NA activity plays a significant role in governing the evolutionary potential of HA in the currently circulating pdmH1N1 lineage.

Evidence of antigenic drift in the fusion machinery core of SARS-CoV-2 spike.

Antigenic drift of SARS-CoV-2 is typically defined by mutations in the N-terminal domain and receptor binding domain of spike protein. In contrast, whether antigenic drift occurs in the S2 domain remains largely elusive. Here, we perform a deep mutational scanning experiment to identify S2 mutations that affect binding of SARS-CoV-2 spike to three S2 apex public antibodies. Our results indicate that spatially diverse mutations, including D950N and Q954H, which are observed in Delta and Omicron variants, respectively, weaken the binding of spike to these antibodies. Although S2 apex antibodies are known to be nonneutralizing, we show that they confer protection in vivo through Fc-mediated effector functions. Overall, this study indicates that the S2 domain of SARS-CoV-2 spike can undergo antigenic drift, which represents a potential challenge for the development of more universal coronavirus vaccines.

Evolution of human H3N2 influenza virus receptor specificity has substantially expanded the receptor-binding domain site.

Hemagglutinins (HAs) from human influenza viruses descend from avian progenitors that bind α2-3-linked sialosides and must adapt to glycans with α2-6-linked sialic acids on human airway cells to transmit within the human population. Since their introduction during the 1968 pandemic, H3N2 viruses have evolved over the past five decades to preferentially recognize human α2-6-sialoside receptors that are elongated through addition of poly-LacNAc. We show that more recent H3N2 viruses now make increasingly complex interactions with elongated receptors while continuously selecting for strains maintaining this phenotype. This change in receptor engagement is accompanied by an extension of the traditional receptor-binding site to include residues in key antigenic sites on the surface of HA trimers. These results help explain the propensity for selection of antigenic variants, leading to vaccine mismatching, when H3N2 viruses are propagated in chicken eggs or cells that do not contain such receptors.

New horizons for comparative studies and meta-analyses.

Comparative analyses and meta-analyses are key tools to elucidate broad biological principles, yet the two approaches often appear different in purpose. We propose an integrated approach that can generate deeper insights into ecoevolutionary processes. Marrying comparative and meta-analytic approaches will allow for (i) a more accurate investigation of drivers of biological variation, (ii) a greater ability to account for sources of non-independence in experimental data, (iii) more effective control of publication bias, and (iv) improved transparency and reproducibility. Stronger integration of meta-analytic and comparative studies can also broaden the scope from species-centric investigations to community-level responses and function-valued traits (e.g., reaction norms). We illuminate commonalities, differences, and the transformative potential of combining these methodologies for advancing ecology and evolutionary biology.

Targeting neuraminidase: the next frontier for broadly protective influenza vaccines.

Current seasonal influenza vaccines, which mainly target hemagglutinin (HA), require annual updates due to the continuous antigenic drift of the influenza virus. Developing an influenza vaccine with increased breadth of protection will have significant public health benefits. The recent discovery of broadly protective antibodies to neuraminidase (NA) has provided important insights into developing a universal influenza vaccine, either by improving seasonal influenza vaccines or designing novel immunogens. However, further in-depth molecular characterizations of NA antibody responses are warranted to fully leverage broadly protective NA antibodies for influenza vaccine designs. Overall, we posit that focusing on NA for influenza vaccine development is synergistic with existing efforts targeting HA, and may represent a cost-effective approach to generating a broadly protective influenza vaccine.

Functional and antigenic characterization of SARS-CoV-2 spike fusion peptide by deep mutational scanning.

The fusion peptide of SARS-CoV-2 spike protein is functionally important for membrane fusion during virus entry and is part of a broadly neutralizing epitope. However, sequence determinants at the fusion peptide and its adjacent regions for pathogenicity and antigenicity remain elusive. In this study, we performed a series of deep mutational scanning (DMS) experiments on an S2 region spanning the fusion peptide of authentic SARS-CoV-2 in different cell lines and in the presence of broadly neutralizing antibodies. We identified mutations at residue 813 of the spike protein that reduced TMPRSS2-mediated entry with decreased virulence. In addition, we showed that an F823Y mutation, present in bat betacoronavirus HKU9 spike protein, confers resistance to broadly neutralizing antibodies. Our findings provide mechanistic insights into SARS-CoV-2 pathogenicity and also highlight a potential challenge in developing broadly protective S2-based coronavirus vaccines.

Stringent and complex sequence constraints of an IGHV1-69 broadly neutralizing antibody to influenza HA stem.

IGHV1-69 is frequently utilized by broadly neutralizing influenza antibodies to the hemagglutinin (HA) stem. These IGHV1-69 HA stem antibodies have diverse complementarity-determining region (CDR) H3 sequences. Besides, their light chains have minimal to no contact with the epitope. Consequently, sequence determinants that confer IGHV1-69 antibodies with HA stem specificity remain largely elusive. Using high-throughput experiments, this study reveals the importance of light-chain sequence for the IGHV1-69 HA stem antibody CR9114, which is the broadest influenza antibody known to date. Moreover, we demonstrate that the CDR H3 sequences from many other IGHV1-69 antibodies, including those to the HA stem, are incompatible with CR9114. Along with mutagenesis and structural analysis, our results indicate that light-chain and CDR H3 sequences coordinately determine the HA stem specificity of IGHV1-69 antibodies. Overall, this work provides molecular insights into broadly neutralizing antibody responses to influenza virus, which have important implications for universal influenza vaccine development.

Trends in Characteristics of Prescription Opioid-related Poisonings among Older Adults in the United States, 2015-2021.

OBJECTIVES: Few studies have considered how trends in opioid poisonings have changed among older adults. The objective of this study was to examine trends in fatal and nonfatal opioid-related poisonings ("exposures") among older adults. METHODS: National poison center data were used to examine trends in characteristics of reported exposures to commonly prescribed opioids between 2015 and 2021 among adults 60 years or older. We estimated the proportion of opioid exposures by demographic characteristics, the specific opioid(s) involved, exposure type, route of administration, other substances co-used, and medical outcomes for each calendar year. We estimated whether there were linear changes in prevalence by year using logistic regression. RESULTS: Although there was a decrease in the number of opioid exposures within the study population from 7706 in 2015 to 7337 in 2021 (a 4.8% decrease, P = 0.04), exposures increased for adults aged 70 to 79 years (a 14.0% increase, P < 0.001). The proportion classified as "abuse" increased by 63.3% ( P < 0.001). There were significant decreases in the proportion involving hydromorphone (a 23.3% decrease, P < 0.001) and morphine (a 22.0% decrease, P < 0.001), with an increase involving buprenorphine (a 216.0% increase, P < 0.001). The proportion increased for co-use of cocaine (a 488.9% increase, P < 0.001) and methamphetamine (a 220.0% increase, P = 0.02), with a decrease in co-use of benzodiazepines (a 25.5% decrease, P < 0.001). The proportion of major medical outcomes increased by 93.9% ( P < 0.001). CONCLUSIONS: National patterns of opioid-related poisonings are shifting among older adults, including the types of opioids involved and co-use of other drugs. These results can inform prevention and harm reduction efforts aimed at older adults.

Leveraging vaccination-induced protective antibodies to define conserved epitopes on influenza N2 neuraminidase.

There is growing appreciation for neuraminidase (NA) as an influenza vaccine target; however, its antigenicity remains poorly characterized. In this study, we isolated three broadly reactive N2 antibodies from the plasmablasts of a single vaccinee, including one that cross-reacts with NAs from seasonal H3N2 strains spanning five decades. Although these three antibodies have diverse germline usages, they recognize similar epitopes that are distant from the NA active site and instead involve the highly conserved underside of NA head domain. We also showed that all three antibodies confer prophylactic and therapeutic protection in vivo, due to both Fc effector functions and NA inhibition through steric hindrance. Additionally, the contribution of Fc effector functions to protection in vivo inversely correlates with viral growth inhibition activity in vitro. Overall, our findings advance the understanding of NA antibody response and provide important insights into the development of a broadly protective influenza vaccine.

Contrasting roles of MERS-CoV and SARS-CoV-2 internal proteins in pathogenesis in mice.

Betacoronaviruses encode an internal (I) gene via an alternative reading frame within the nucleocapsid gene, called ORF8b for Middle-East respiratory syndrome coronavirus (MERS-CoV) and ORF9b for severe acute respiratory syndrome coronavirus (SARS-CoV) and SARS-CoV-2. Previous reports suggested that proteins 8b and 9b are involved in evading multiple innate immune signaling pathways. However, their roles in mediating pathogenesis in infected animals have not been determined. In this study, we abrogated the expression of protein 8b in MERS-CoV and protein 9b in SARS-CoV-2. Using mouse models of MERS-CoV and SARS-CoV-2 infection, we found that MERS-CoV lacking protein 8b expression was more virulent, while SARS-CoV-2 lacking protein 9b expression was attenuated compared with the respective wild-type viruses. Upon further analysis, we detected increased levels of type I interferon and enhanced infiltration of immune cells to the lungs of mice infected with MERS-CoV lacking protein 8b expression. These data suggest that the I protein of MERS-CoV plays a role in limiting pathogenesis while that of SARS-CoV-2 enhances disease severity. IMPORTANCE The function of betacoronavirus internal protein has been relatively understudied. The earliest report on the internal protein of mouse hepatitis virus suggested that the internal protein is a structural protein without significant functions in virus replication and virulence. However, the internal proteins of severe acute respiratory syndrome coronavirus (SARS-CoV), Middle-East respiratory syndrome coronavirus, and SARS-CoV-2 have been shown to evade immune responses. Despite the reported functions of the internal protein in these highly pathogenic human coronaviruses, its role in mediating pathogenesis in experimentally infected animals has not been characterized. Our data indicated that despite the similar genomic location and expression strategy of these internal proteins, their effects on virulence are vastly different and virus specific, highlighting the complexity between host-virus interaction and disease outcome.

Widespread impact of immunoglobulin V-gene allelic polymorphisms on antibody reactivity.

The ability of the human immune system to generate antibodies to any given antigen can be strongly influenced by immunoglobulin V-gene allelic polymorphisms. However, previous studies have provided only limited examples. Therefore, the prevalence of this phenomenon has been unclear. By analyzing >1,000 publicly available antibody-antigen structures, we show that many V-gene allelic polymorphisms in antibody paratopes are determinants for antibody binding activity. Biolayer interferometry experiments further demonstrate that paratope allelic polymorphisms on both heavy and light chains often abolish antibody binding. We also illustrate the importance of minor V-gene allelic polymorphisms with low frequency in several broadly neutralizing antibodies to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and influenza virus. Overall, this study not only highlights the pervasive impact of V-gene allelic polymorphisms on antibody binding but also provides mechanistic insights into the variability of antibody repertoires across individuals, which in turn have important implications for vaccine development and antibody discovery.

An explainable language model for antibody specificity prediction using curated influenza hemagglutinin antibodies.

Despite decades of antibody research, it remains challenging to predict the specificity of an antibody solely based on its sequence. Two major obstacles are the lack of appropriate models and inaccessibility of datasets for model training. In this study, we curated a dataset of >5,000 influenza hemagglutinin (HA) antibodies by mining research publications and patents, which revealed many distinct sequence features between antibodies to HA head and stem domains. We then leveraged this dataset to develop a lightweight memory B cell language model (mBLM) for sequence-based antibody specificity prediction. Model explainability analysis showed that mBLM captured key sequence motifs of HA stem antibodies. Additionally, by applying mBLM to HA antibodies with unknown epitopes, we discovered and experimentally validated many HA stem antibodies. Overall, this study not only advances our molecular understanding of antibody response to influenza virus, but also provides an invaluable resource for applying deep learning to antibody research.

Stringent and complex sequence constraints of an IGHV1-69 broadly neutralizing antibody to influenza HA stem.

IGHV1-69 is frequently utilized by broadly neutralizing influenza antibodies to the hemagglutinin (HA) stem. These IGHV1-69 HA stem antibodies have diverse complementarity-determining region (CDR) H3 sequences. Besides, their light chains have minimal to no contact with the epitope. Consequently, sequence determinants that confer IGHV1-69 antibodies with HA stem specificity remain largely elusive. Using high-throughput experiments, this study revealed the importance of light chain sequence for the IGHV1-69 HA stem antibody CR9114, which is the broadest influenza antibody known to date. Moreover, we demonstrated that the CDR H3 sequences from many other IGHV1-69 antibodies, including those to HA stem, were incompatible with CR9114. Along with mutagenesis and structural analysis, our results indicate that light chain and CDR H3 sequences coordinately determine the HA stem specificity of IGHV1-69 antibodies. Overall, this work provides molecular insights into broadly neutralizing antibody responses to influenza virus, which have important implications for universal influenza vaccine development.

Defining and Studying B Cell Receptor and TCR Interactions.

BCRs (Abs) and TCRs (or adaptive immune receptors [AIRs]) are the means by which the adaptive immune system recognizes foreign and self-antigens, playing an integral part in host defense, as well as the emergence of autoimmunity. Importantly, the interaction between AIRs and their cognate Ags defies a simple key-in-lock paradigm and is instead a complex many-to-many mapping between an individual's massively diverse AIR repertoire, and a similarly diverse antigenic space. Understanding how adaptive immunity balances specificity with epitopic coverage is a key challenge for the field, and terms such as broad specificity, cross-reactivity, and polyreactivity remain ill-defined and are used inconsistently. In this Immunology Notes and Resources article, a group of experimental, structural, and computational immunologists define commonly used terms associated with AIR binding, describe methodologies to study these binding modes, as well as highlight the implications of these different binding modes for therapeutic design.

Widespread impact of immunoglobulin V gene allelic polymorphisms on antibody reactivity.

The ability of human immune system to generate antibodies to any given antigen can be strongly influenced by immunoglobulin V gene (IGV) allelic polymorphisms. However, previous studies have provided only a limited number of examples. Therefore, the prevalence of this phenomenon has been unclear. By analyzing >1,000 publicly available antibody-antigen structures, we show that many IGV allelic polymorphisms in antibody paratopes are determinants for antibody binding activity. Biolayer interferometry experiment further demonstrates that paratope allelic mutations on both heavy and light chain often abolish antibody binding. We also illustrate the importance of minor IGV allelic variants with low frequency in several broadly neutralizing antibodies to SARS-CoV-2 and influenza virus. Overall, this study not only highlights the pervasive impact of IGV allelic polymorphisms on antibody binding, but also provides mechanistic insights into the variability of antibody repertoires across individuals, which in turn have important implications for vaccine development and antibody discovery.

Sublethal consequences of ultraviolet radiation exposure on vertebrates: Synthesis through meta-analysis.

Ultraviolet radiation (UVR) from the sun is a natural daytime stressor for vertebrates in both terrestrial and aquatic ecosystems. UVR effects on the physiology of vertebrates manifest at the cellular level, but have bottom-up effects at the tissue level and on whole-animal performance and behaviours. Climate change and habitat loss (i.e. loss of shelter from UVR) could interact with and exacerbate the genotoxic and cytotoxic impacts of UVR on vertebrates. Therefore, it is important to understand the range and magnitude of effects that UVR can have on a diversity of physiological metrics, and how these may be shaped by taxa, life stage or geographical range in the major vertebrate groups. Using a meta-analytical approach, we used 895 observations from 47 different vertebrate species (fish, amphibian, reptile and bird), and 51 physiological metrics (i.e. cellular, tissue and whole-animal metrics), across 73 independent studies, to elucidate the general patterns of UVR effects on vertebrate physiology. We found that while UVR's impacts on vertebrates are generally negative, fish and amphibians were the most susceptible taxa, adult and larvae were the most susceptible life stages, and animals inhabiting temperate and tropical latitudes were the most susceptible to UVR stress. This information is critical to further our understanding of the adaptive capacity of vulnerable taxon to UVR stress, and the wide-spread sublethal physiological effects of UVR on vertebrates, such as DNA damage and cellular stress, which may translate up to impaired growth and locomotor performance. These impairments to individual fitness highlighted by our study may potentially cause disruptions at the ecosystem scale, especially if the effects of this pervasive diurnal stressor are exacerbated by climate change and reduced refuge due to habitat loss and degradation. Therefore, conservation of habitats that provide refuge to UVR stress will be critical to mitigate stress from this pervasive daytime stressor.

High-throughput identification of prefusion-stabilizing mutations in SARS-CoV-2 spike.

Designing prefusion-stabilized SARS-CoV-2 spike is critical for the effectiveness of COVID-19 vaccines. All COVID-19 vaccines in the US encode spike with K986P/V987P mutations to stabilize its prefusion conformation. However, contemporary methods on engineering prefusion-stabilized spike immunogens involve tedious experimental work and heavily rely on structural information. Here, we establish a systematic and unbiased method of identifying mutations that concomitantly improve expression and stabilize the prefusion conformation of the SARS-CoV-2 spike. Our method integrates a fluorescence-based fusion assay, mammalian cell display technology, and deep mutational scanning. As a proof-of-concept, we apply this method to a region in the S2 domain that includes the first heptad repeat and central helix. Our results reveal that besides K986P and V987P, several mutations simultaneously improve expression and significantly lower the fusogenicity of the spike. As prefusion stabilization is a common challenge for viral immunogen design, this work will help accelerate vaccine development against different viruses.

Imprinted Anti-Hemagglutinin and Anti-Neuraminidase Antibody Responses after Childhood Infections of A(H1N1) and A(H1N1)pdm09 Influenza Viruses.

Immune imprinting is a driver known to shape the anti-hemagglutinin (HA) antibody landscape of individuals born within the same birth cohort. With the HA and neuraminidase (NA) proteins evolving at different rates under immune selection pressures, anti-HA and anti-NA antibody responses since childhood influenza virus infections have not been evaluated in parallel at the individual level. This is partly due to the limited knowledge of changes in NA antigenicity, as seasonal influenza vaccines have focused on generating neutralizing anti-HA antibodies against HA antigenic variants. Here, we systematically characterized the NA antigenic variants of seasonal A(H1N1) viruses from 1977 to 1991 and completed the antigenic profile of N1 NAs from 1977 to 2015. We identified that NA proteins of A/USSR/90/77, A/Singapore/06/86, and A/Texas/36/91 were antigenically distinct and mapped N386K as a key determinant of the NA antigenic change from A/USSR/90/77 to A/Singapore/06/86. With comprehensive panels of HA and NA antigenic variants of A(H1N1) and A(H1N1)pdm09 viruses, we determined hemagglutinin inhibition (HI) and neuraminidase inhibition (NI) antibodies from 130 subjects born between 1950 and 2015. Age-dependent imprinting was observed for both anti-HA and anti-NA antibodies, with the peak HI and NI titers predominantly detected from subjects at 4 to 12 years old during the year of initial virus isolation, except the age-independent anti-HA antibody response against A(H1N1)pdm09 viruses. More participants possessed antibodies that reacted to multiple antigenically distinct NA proteins than those with antibodies that reacted to multiple antigenically distinct HA proteins. Our results support the need to include NA proteins in seasonal influenza vaccine preparations. IMPORTANCE Seasonal influenza vaccines have aimed to generate neutralizing anti-HA antibodies for protection since licensure. More recently, anti-NA antibodies have been established as an additional correlate of protection. While HA and NA antigenic changes occurred discordantly, the anti-HA and anti-NA antibody profiles have rarely been analyzed in parallel at the individual level, due to the limited knowledge on NA antigenic changes. By characterizing NA antigenic changes of A(H1N1) viruses, we determined the anti-HA and anti-NA antibody landscape against antigenically distinct A(H1N1) and A(H1N1)pdm09 viruses using sera of 130 subjects born between 1950 and 2015. We observed age-dependent imprinting of both anti-HA and anti-NA antibodies against strains circulated during the first decade of life. A total of 67.7% (88/130) and 90% (117/130) of participants developed cross-reactive antibodies to multiple HA and NA antigens at titers ≥1:40. With slower NA antigenic changes and cross-reactive anti-NA antibody responses, including NA protein in influenza vaccine preparation may enhance vaccine efficacy.

Mutational fitness landscape of human influenza H3N2 neuraminidase.

Influenza neuraminidase (NA) has received increasing attention as an effective vaccine target. However, its mutational tolerance is not well characterized. Here, the fitness effects of >6,000 mutations in human H3N2 NA are probed using deep mutational scanning. Our result shows that while its antigenic regions have high mutational tolerance, there are solvent-exposed regions with low mutational tolerance. We also find that protein stability is a major determinant of NA mutational fitness. The deep mutational scanning result correlates well with mutational fitness inferred from natural sequences using a protein language model, substantiating the relevance of our findings to the natural evolution of circulating strains. Additional analysis further suggests that human H3N2 NA is far from running out of mutations despite already evolving for >50 years. Overall, this study advances our understanding of the evolutionary potential of NA and the underlying biophysical constraints, which in turn provide insights into NA-based vaccine design.